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Problem  Statement: 

Matrix  eigenvalue  |)i'ol)lems  play  a  significani  roh'  in  dynamic  analysis  of  structures.  The  natural 
fre([uencies  and  the  modes  oi  a.  free  vibrating  system  are  intimately  related  to  the  eigenvMues 
and  the  eigenvectors  of  tlie  geneialized  system  of  ( haracteristic  ecpiations  correspondiirg  to  that 
vibrating  svstem.  To  calculate  even  a  few  ol  the  natural  frecpiencies  and  modes  of  a  large  scale 
structure  is  computationally  intensive,  k  frecpiently  encountered  problem  in  structural  dynamics 
is  how  to  take  into  account  analysis  and  design  changes  introduced  after  the  structural  dynamic 
analysis  has  been  completed  and  the  natural  fre(|uenries  and  modes  have  been  computed.  The 
problem  is  how  the  natural  frequencies  and  modes  obtained  from  a  previous  analysis  can  be  used  to 
derive  the  response  of  the  new  modified  structure,  without  extensive  additional  computations.  In 
particular,  we  focus  on  tlie  structural  dynamic  re-analysis  prol)lem  that  the  system  is  modified  in  a 
few  small  localized  regions  the  incremental  stiffness  matrix  and  the  incremental  mass  matrix  are 
low  rank  matrices,  but  t  heir  magnit  ude  need  not  In*  small  relative  to  the  original  matrices.  Dynamic 
solution  procedures  that  are  suitable  for  parallel  computer  environments  are  also  investigated. 

Summary  of  Result? 

In  this  project,  we  hav(’  developed  and  impletiiented  a  number  of  parallel  algorithms  related  to 
matrix  factorization  and  generalized  eigenvalue  problems.  Furthermore,  a  Lanczos  based  method 
for  the  modified  eigenvalue  problem  has  been  develojred.  The  following  summarizes  the  results  of 
this  project: 

1.  Development  and  implementation  of  a  new  |)arallel  sparse  factorization  method(Ref.  1,3). 

Specific  contribution  includes  au  efficient  procedure  for  the  solution  of  triangular  systems, 
which  is  generallv  known  as  one  of  the  most  difficult  problems  for  the  parallelization  of  direct 
solution  methods.  The  parallel  j)rocedure  for  the  solution  of  triangular  systems  is  particularly 
important  for  generalized  eigeinnlue  problem  which  requires  solution  of  linear  systems  at  each 
Lanczos  iteration  stej).  □ 

2.  A  parallel  implementation  of  a  robust  Lanczos  algorithm  for  generalized  eigenvalue  problems  . 

(Ref.  1,10).  The  parallel  algorithm  has  been  implemented  in  a  finite  element  structural  anal- 

ysis  program.  To  ensure  the  lobustness  and  accuracy,  the  parallel  procedure  implemented  - 

includes  spectral  shifts,  partial  and  selective  re-orthogonalizations.  We  have  also  tested  vari- 


ous  approaches 

and  develo|K'd  a  robust  implementation  so  that  the  finite  element  matrices  can 

be  in  the  form 

of  ('ither  element  matrices  or  assembled  system  matrices.  Tradeoffs  between 

_ _  j 

spectral  shifts. 

factorization  and  re-orthogonaliza.tion  have  been  studied. 

.} 
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The  development  of  a  now  method  for  solving  perturbed  symmetric  eigenvalue  problems 
(Ref.  3.5.9).  This  method  is  Intsed  on  an  extension  of  the  block  Lanezos  method  for  standard 
and  generalized  eigt'iivalne  problems.  Ihis  lu'w  method  is  capable  ol  updating  the  modified 
eigenvalues  and  eigenvectors  e\'en  when  the  modification  made  to  the  structures  may  cause 
significant  changes  to  the  mode  shapes  ol  the  system. 

In  addition,  we  have  also  studied  the  following  problems  related  to  eigenvalue  problems  and  the 
Lanezos  and  conjugate  gradient  methods  lor  structural  analysis; 

1.  A  study  of  the  solution  methods  for  the  minimization  of  the  generalized  Rayleigh  Quotient 
subject  to  linear  and  cpiadratic  constraints  (Ref.  7).  Applications  of  the  methods  for  the 
evaluation  of  thin  shell  finite  elements  are  illustrated. 

2.  Development  of  a  robust  and  stable  incomplete  (  holesky  (IC  )  lactorization  (Rel.  2,8).  In 
particular,  we  ha\('  examined  in  detail  the  trath'offs  between  the  storage  reejnirement  and  the 
effectiveness  (in  tennis  of  the  total  CPL  time  for  the  solution)  of  the  IC.  pieconditioneis  foi 
the  conjugate  gradient  method.  .\  techniepie  has  been  developed  to  manage  both  the  storage 
space  and  to  ()|rtimiz('  the  "goodness’  ol  tin'  pre'conditioning  matrix. 

3.  Development  of  a  simple  bni  elfi'ctive  procedure  to  enhance  the  Irlock  conjugate  gradient 
method  for  the  solution  of  a  linear  system  of  equations  with  multiple  right  hand  sides  (Ref. 
2,6).  This  development  is  important  since  engineering  structures  are  often  designed  to  resist 
multi  [lie  load  cases. 

4.  A  study  of  the  use  of  the  conjugate  gradient  method  for  nonlinear  analysis  of  structures 
(including  snap  through  and  simple  bilurcation  ol  structures)  (Ref.  2).  Based  on  the  rela¬ 
tionship  between  the  conjugate  gradient  and  the  Lanezos  methods,  we  have  examined  how  to 
take  advantage  of  tlie  information  (such  as  the  eigenvalues,  the  residual  and  search  direction 
vectors)  in  soK'ing  nonlinear  problems.  For  example,  by  simply  keeping  track  of  the  inertia 
of  the  (tangential)  stiffness  matrix,  we  can  determine  easily  the  secondary  load-displacement 
path  related  to  simiile  notdinear  bifurcation  ol  structures. 

5.  A  comparison  study  on  a  variety  of  preconditioners  (including  polynomial  preconditioners, 
incomplete  Cholesky  factors,  domain  decomposition  and  a  combination  of  the  preconditioning 
methods)  for  conjugate  gradient  methods  on  a  distributed  memory  parallel  computer  (Ref. 
2).  Both  linear  and  nonlinear  structural  problems  have  been  exanrined  in  this  study. 

All  parallel  solution  procedures  have  been  implemented  and  tested  on  an  Intel  iS60  parallel  com¬ 
puter. 
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